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PERFORMANCE  OF  CELLULAR  FREQUENCY-HOPPED 
SPREAD-SPECTRUM  RADIO  NETWORKS 


1.  INTRODUCTION 

In  many  problems  of  practical  interest  in  communication  networks  it  is  advantageous 
to  combine  the  use  of  spread-spectrum  signaling  techniques  with  networking  techniques 
in  order  to  provide  enhanced  multi-user  capabibties  and  a  higher  degree  of  resistance  to 
interference,  whether  hostile  or  benign. 

However,  many  of  the  spread-spectrum  network  models  proposed  in  the  past  have 
lacked  precision  in  their  descriptions  of  the  effects  of  spread-spectrum  techniques  on  net¬ 
work  performance.  The  general  practice  in  papers  like  [1]  and  [3]  has  been  to  model  the 
effects  of  spread-spectrum  techniques  via  the  uf*>  of  some  threshold  (a  number  of  users 
beyond  which  communication  is  not  possible)  or  via  some  processing  gain.  Unfortunately, 
such  approaches  generally  do  not  provide  accurate  models,  and  the  results  obtained  may  be 
optimistic  or  pessimistic,  depending  upon  assumptions  made.  In  addition,  these  analyses 
do  not  accurately  account  for  data  modulation  and  forward  error-control  coding. 

In  [1],  a  model  was  presented  for  a  frequency-hopped  (FH)  spread-spectrum  multi¬ 
ple  access  (SSMA)  digital  cellular  telephone  network  and  an  expression  was  derived  for 
throughput.  Notable  features  of  that  model  were  the  characterization  of  user  mobility 
in  terms  of  a  two-dimensional  Poisson  process  and  the  treatment  of  transmitter  power 
attenuation  with  distance.  This  model  superceded  that  of  [18]. 

In  [2],  the  authors  described  and  analyzed  a  digital  cellular  mobile  network  using  some 
of  the  techniques  of  [1].  While  adding  to  the  model  the  interference  due  to  users  in  other 

ce"'  the  authors  L-ft  out  the  spatial  Poisson  model  and  did  not  account  explicitly  for 
Mar.uacrip*  anorov-d  A.ig-^r,;  3u.  i9d9. 
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power  attentuation  with  distance.  Spatial  Poisson  modeling  was  also  used  in  [3],  in  which 
the  authors  modeled  spread-spectrum  multihop  networks.  In  [4],  a  more  precise  approach 
was  taken  to  the  spread-spectrum  multihop  network  problem,  while  using  some  similar 
network  features  in  the  analysis. 

In  this  paper,  the  modeling  techniques  described  above  are  combined  with  a  precise 
characterization  of  the  effects  of  frequency  hopping  for  the  purpose  of  modeling  frequency- 
hopped  digital  cellular  networks.  This  analysis  also  incorporates  different  data  modulation 
types  (coherent  and  noncoherent)  and  forward  error-control  coding.  Also,  unlike  several 
of  the  previous  analyses,  some  of  the  techniques  presented  here  allow  for  the  computation 
of  bit/symbol  or  packet  error  probability,  as  desired.  In  addition  to  modeling  a  ceU  of  a 
FH/SSMA  cellular  mobile  radio  network,  this  also  corresponds  to  a  more  accurate  analysis 
of  communications  in  the  presence  of  secondary  multiple-access  interference  in  that  it 
accounts  for  varying  powers  for  different  interferers  according  to  their  distances  from  the 
receiver  (compare,  for  example,  to  [11]  and  [19]). 

An  overview  of  the  contents  of  this  paper  is  as  follows.  In  Section  II,  the  network 
model  is  presented.  Section  III  contains  analyses  of  the  model  for  coherent  and  non¬ 
coherent  modulation/demodulation  using  approximation  techniques.  Section  IV  provides 
an  exact  analysis  for  coherent  systems  using  a  characteristic  function  method.  In  Section 
V,  numerical  results  are  presented  and  compared  with  previously  derived  results.  Section 
6  contains  the  concluding  remarks. 

2.  NETWORK  MODEL 

The  network  model  used  in  this  paper  is  similar  to  the  one  found  in  [1].  Users  transmit 
packets  to  other  users  via  fixed  relay  (or  “base”)  stations.  As  is  the  case  with  [1],  we  will 
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focus  on  the  user-to-station  channel.  Analysis  of  the  station-to-user  channel  would  then 
proceed  in  a  similar  fashion 

Let  the  number  of  users  per  unit  area  be  Poisson  distributed  with  parameter  A.  Also, 
let  M  =  {x,y)  represent  a  position  on  the  plane.  Denote  the  origin  by  0  =  (0,0);  this  will 
be  assumed,  without  loss  of  generality,  to  be  the  position  of  the  fixed  station. 

Systems  like  this  generally  offer  voice  and  data  service.  In  the  case  of  voice  service,  a 
steady  stream  of  information  must  be  transmitted  and  delay  can  not  be  tolerated;  however, 
moderate  error  rates  can  be  tolerated.  In  the  case  of  data  service,  traffic  is  bursty  and  large 
delays  can  be  tolerated  but  small  error  rates  are  necessary.  In  this  paper,  it  is  assumed 
that  voice  (telephone)  service  is  predominant. 

Transmitter-oriented  assignment  of  hopping  patterns  is  also  assumed.  Each  mobile 
has  its  own  frequency-hopping  pattern  for  transmission  and  reception;  the  bcise  station 
has  a  list  of  all  of  the  patterns  and  can  listen  to  several  of  them  (or  even  aU  of  them) 
simultaneously.  We  restrict  our  attention  to  the  case  of  mobile-to-base  communication. 
We  consider  the  situation  in  which  a  mobile  station  has  already  established  communication 
with  the  base  station  at  the  center  of  the  cell  and  now  communicates  in  the  presence  of 
secondary  (or  multiple-access)  interference  from  other  mobiles.  It  is  assumed  that  only 
the  other  terminals  in  the  cell  of  interest  may  interfere,  which  reflects  the  use  of  other 
frequency  bands  in  all  adjacent  cells,  and  that  frequencies  are  only  reused  in  cells  distant 
enough  so  that  the  effects  are  negligible. 

Synchronization  at  the  packet  level  is  assumed  feasible  for  all  users.  Thus,  uncertainties 
in  the  timing  between  different  users  are  small  relative  to  packet  duration;  however,  since 
they  might  not  be  small  relative  to  the  dwell  time  of  the  frequency  hopper,  the  system  is 
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modeled  as  a  general  asynchronous  frequency  hopping  system. 

3.  APPROXIMATIONS  AND  BOUNDS  ON  PACKET  ERROR 
PROBABILITY  AND  THROUGHPUT 

Our  approximation  methods  are  developed  below.  These  can  be  used  for  coherently 
and  noncoherently  modulated  systems  with  or  without  forward  error-control  coding. 

3.1.  Analysis  Based  on  Gaussian  Approximation  Techniques 

3.1.1.  Noncoherent  Signaling 

We  begin  with  noncoherent  signaling.  In  this  case,  the  system  is  assumed  to  utilize  M- 
ary  frequency-shift  keying  (MFSK);  binary  FSK  is  a  special  case  of  MFSK,  and  therefore 
the  analysis  holds  for  this,  as  well.  Noncoherent  demodulation  is  used  at  the  receiver. 

The  signal  transmitted  by  the  k-th  mobile  user  is  of  the  form 

Skit)=  cos{2ir[f,  +  bk{t)A  + fk{t)]t  +  0kit)}  (1) 

where  Pq  is  the  transmitted  signal  power  and  is  assumed  to  be  the  same  for  all  mobile 
stations;  bk{t)  =  where  e  {1,2, ...,M}  is  the  A:-th  sequence  of  M-ary  infor¬ 

mation  symbols  of  duration  T  each;  ’^(t)  is  a  pulse-shaping  waveform;  0k{t)  is  the  sum  of 
phase  shifts  due  to  modulation  and  frequency  hopping  generated  by  the  k-th  transmitter; 
fc  is  the  carrier  frequency;  2A  is  the  spacing  between  any  two  consecutive  tones  out  of 
the  M  different  tones;  and  fk{t)  is  the  k-th  hopping  sequence. 

In  our  analysis,  we  consider  only  the  attenuation  of  signal  power  with  distance.  We 
assume  that  the  attenuation  and  the  noise  spectral  density  are  constant  over  the  hopping 
bandwidth.  In  many  practical  systems  attenuation  varies  over  the  hopping  bandwidth. 
It  would  be  possible  to  include  such  frequency-dependent  effects  in  our  model;  however. 


to  do  so  would  further  complicate  an  already  complex  and  computationally  demanding 
calculation.  * 

Let  P(M)  denote  the  received  power  of  a  signal  originating  at  the  mobile  station 
located  at  the  point  M  of  the  plane  and  received  at  the  origin.  We  will  assume  that  P{M) 
is  given  by 


P(M)  =  { 


if  |0M|  <  ro 


(2) 


[  Po|0Ml-“,  if  lOM!  >  ro 

where  |0M1  denotes  the  Euclidean  distance  between  M  and  the  origin  0,  ro  is  some  small 
radius  within  which  P{M)  is  assumed  to  be  approximately  constant,  Pq  is  the  transmitted 
power,  and  q  is  a  constant  that  reflects  the  speed  of  attenuation,  ro  is  necessary  because, 
as  M  is  brought  arbitrary  close  to  0,  |0M1"“  approaches  infinity.  A  particular  signal 
arriving  at  the  receiver  (and  originating  from  Mfe)  ha.s  the  form 


rfc(t)  =  ^P{Mk)^{t)  cos  [f,  -f  6it(0A  +  A(t)]  t  +  (3) 

where  sj2P(Mk)  represents  the  received  signal  amplitude  and  OkW  is  the  total  pha^e  shift. 
Since  the  system  under  consideration  is  a  multiple-access  one,  the  received  signal  has  the 
form 

^(0  =  Yi  -  Tk)  +  nit)  (4) 

k 

where  each  Vkit)  is  given  by  (3),  Tk  is  a  random  delay,  and  n{t)  is  zero-mean  additive  white 
Gaussian  noise  (AWGN)  with  two-sided  spectral  density  Ao/2. 

'It  would  actually  involve  deriving  the  bit  or  packet  error  probability  (accounting  for  both  power 
and  frequency  losses  and  other-user  interference)  conditioned  on  the  system  operating  at  a  particular 
hopping  frequency  and  then  integrating  over  the  probability  density  function  of  the  system  visiting  each 
hopping  frequency.  Therefore,  one  additional  integration  would  be  necessary  besides  the  integrations 
already  performed  in  our  subsequent  andysis. 
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Fig.  1  shows  the  receiver  for  the  tth  signal  in  a  binary  FSK  system;  the  A/-ary  receiver 
is  similar  in  form  but  w'ith  M  branches.  The  ith  signal  is  the  signal  of  interest.  r(t)  is 
dehopped  and  the  dehopper  output  is  given  by 

X  cos  [27r[/c +  &*(<- (5) 

Here  /*(•)  and  /,(■)  represent  the  fcth  and  ith  hopping  sequences,  respectively;  nj{t)  is  a 
zero-mean  Gaussian  random  process  with  spectral  density  A'o/8;  ^(u,r)  =  1,  if  u  =  r, 
and  0,  otherwise;  and  (t>k{t)  represents  the  total  phase  shift  of  the  signal  after  dehopping. 
The  sufficient  statistics  are  formed  as  shown  in  Fig.  1.  Were  it  not  for  multiple- 
access  interference,  the  computation  of  the  probability  of  error  for  the  system  described 
above  (a  binary  FSK  system  with  noncoherent  demodulation)  would  be  straightforward, 
as  for  example  in  [10].  The  presence  of  other-user  interference  is  due  to  the  partial  overlap 
of  the  frequency-hopping  patterns  assigned  to  the  different  mobiles;  this  is  termed  sec¬ 
ondary  multiple-access  interference.  By  contrast,  the  interference  generated  when  a  num¬ 
ber  of  transmitters  using  the  same  hopping  pattern  attempt  to  access  the  same  receiver  is 
termed  primary  interference.  To  incorporate  the  effects  of  secondary  other-user  interfer¬ 
ence  is  difficult  because  it  brings  into  the  computation  of  the  error  probability  features  of 
each  transmitting  signal  (i.e.,  time  delays,  phase  angles,  and  data  streams)  which  have  a 
statistical  nature  and  consequently  a  multi-dimensional  integration  needs  to  be  performed 
to  evaluate  the  average  of  the  error  probability  with  respect  to  all  these  variables.  This 
difficulty  can  be  overcome  by  means  of  conditioning  the  error  probability  on  the  potential 
number  of  interfering  users  and  the  actual  number  of  interfering  users  causing  hits  on  the 
desirable  received  signal.  Recall  that  a  hit  occurs  when  fk{t  -  Tk)  =  /,(<)  (refer  to  [5]  and 
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[6]);  if  this  equality  is  satisfied  for  the  entire  duration  of  a  particular  hop  (dwoll-time)  of 
the  i-th  user  we  say  that  user  k  causes  a  full  hit  during  that  hop,  otherwise  we  say  that  the 
Ar-th  user  causes  a  partial  hit  to  the  i-th  user.  To  evaluate  the  conditional  symbol-error 
probability  just  described,  we  use  the  Gaussian  approximation  technique,  as  developed  in 
[5]  and  [6],  where  the  Zc,m  and  Z,,m  representing  the  outputs  of  the  matched  filters  of 
the  in-phase  and  quadrature  parts  of  branch  m  (m  takes  values  -f  1  and  —1)  of  the  BFSK 
demodulator  of  Fig.  1  were  approximated  by  zero-mean  Gaussian  random  variables  with 
the  same  second-order  moments.  This  Gaussiar  approximation  is  one  of  the  most  popular 
tools  for  computing  the  other-user  interference  in  multi-user  spread-spectrum  systems  (see 
[5], [6],  and  [17])  and  has  satisfactory  accuracy  in  most  situations  as  the  comparisons  with 
the  exact  expressions  for  the  error  probability  in  [6]  and  [17]  attest  to.  Given  this,  the 
symbol  error  probability  can  be  derived  (see  [5]  or  [6])  and  is  given  by 


are  full  hits.  Note  that  this  causes  the  result  to  be  pessimistic,  since  we  assume  that  the 
interference  caused  by  any  hit  is  that  of  a  full  hit,  while  the  probability  of  a  hit  is  the  sum 
of  the  probabilities  of  fuU  and  partial  hits.  Under  this  assumption  (7)  becomes 

4  =  hl^&dOA/.Dlog^A/l  '  +  (8) 

in  which  is  used  to  denoie  the  variance  of  interference  under  the  assumption  that  any 
hit  that  occurs  is  a  fuD  hit  (note  that  this  is  independent  of  which  of  the  other  signals 
is  interfering).  Expressions  for  aj  can  be  found  in  [5],  in  whicii  it  is  also  shown  that,  for 
random  FH  patterns,  (7^  —  m^./Af.  For  a  rectangular  shaping  waveform  ^'(<),  m,^,  =  1/3; 
for  a  sine  shaping  waveform,  m^,  =  (15  +  7r^)/127r^. 

At  this  point  we  use  a  three-term  exponential  approximation  to  P^,,,  the  symbol  cor¬ 
rectness  vrobability.  An  exponential  approximation  considerably  facilitates  the  remainder 
of  this  analysis;  this  is  justified  by  the  fact  that  terms  which  are  exponential  in  the  number 
of  interfering  users  can  be  efficiently  averaged  (can  be  actually  put  in  closed  form)  with 
respect  to  Poisson  distributions  on  the  number  of  potentiaJly  interfering  users.  One  could 
use  an  arbitrary  number  of  terms  in  this  exponential  approximation  but  we  found  out 
that  the  complexity  of  the  necessary  subsequent  operations  (for  taking  averages)  would 
be  prohibitive  if  the  number  oi  terms  were  larger  than  three.  Moreover,  since  very  satis¬ 
factory  accuracy  (curve-fitting)  can  be  achieved  with  three  terms,  we  used  three  terms  in 
our  analysis.  Such  an  approximation  can  be  made  using  techniques  like  Prony's  method 
(see  [8]),  or  curve-fitting.  A  least-squares  curve  fit  was  used  to  obtain  the  results  shown 
below;  as  a  result,  we  have 

Pc.,  =  1-Pe..(3-) 
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(9) 


present,  and  on  the  number  of  hits  that  actually  occur  given  the  above  number  of  potential 
interferers.  However,  since  our  ultimate  goal  is  to  determine  packet  throughput,  w-e  require 
a  packet  error  correctness  probability.  This  is  evaluated  in  Section  3.1.3  below. 
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3.1.2.  Coherent  Signaling 

The  development  of  the  probability  expression  that  corresponds  to  (11)  above  for  the 
coherent  case  is  similar  to  the  one  carried  out  above  for  the  noncoherent  case.  In  this 
case,  the  system  is  assumed  to  utilize  binary  phase-shift  keying  (BPSK)  with  coherent 
demodulation.  The  signal  transmitted  by  the  fcth  mobile  station  is  of  the  form 

Skit)  =  \/^bkit)^it)  cos[2n[fc  +  fkit)]t  +  9kit)]  (12) 

in  which  aU  quantities  are  as  above,  except  that  6  {  —  1,  -1-1}  (we  now  have  information 
bits  rather  than  M-ary  symbols).  From  (12)  we  obtain  the  following  expression  for  the 
/:th  received  symbol 

Tkit)  =  yj2P{Mk)bkimit)  cos[27r[/c  +  A(0]<  +  ^fc(0]  •  (13) 

The  received  signal,  given  by  (4),  goes  through  the  same  dehopper  and  the  dehopper 
output  is  given  by 

^<^(0  =  T.^\[^^Hhit-Tk),f,{mk{t-rk)^it-Tk) 

cos[2Trfct  +  <f>kit)]}  +  nd{t).  (14) 


This  is  fed  into  a  coherent  demodulator  and  the  sufficient  statistic  for  the  reception  of  the 

/th  bit  of  the  ith  transmitted  signal  is  given  by 

(;+i)T 

Z,  =  J  rd(i)’i/{t)  cos(27r/cOdt.  (15) 

IT 

Assuming  that  the  receiver  is  synchronized  to  the  tth  signal,  Z,  can  be  rewritten  as  ‘ 


Z.  = 


T  6' 


(») 


+  E 


h,, 


-t-  fi 


(16) 
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where  Ik  ,  represents  the  normalized  interference  caused  by  the  fcth  sign^il  to  the  reception 
of  the  ith  signal  and  h  is  a  zero-mean  Gaussian  random  variable  with  variance  AqT/IG. 
Refer  to  [1]  for  a  detailed  description  of  the  receiver  of  a  coherent  hybrid  direct-sequence 
frequency-hopped  system.  We  only  need  eliminate  the  direct-sequence  part  and  set  A"  =  1 
for  the  number  of  chips  per  bit  to  obtain  the  coherent  purely  frequency-hopped  system. 

The  decision  process  is  concerned  with  determining  whether  the  current  bit  is  -)-l  or 
-1.  The  test  is,  therefore,  to  determine  whether  Z,  >  0  or  Z,  <  0.  This  allows  some 
further  manipulation  of  Z,.  Z,.  the  modified  test  statistic  after  we  divide  by  Pq,  is  of  the 
form 

=  ^J9i\OM,\)bl'^  -h  Yl\/9i\OMk\)Ik,,  +  n  (17) 

where  n  is  a  Gaussian  random  variable  with  variance  No/2Eb,  where  Eb  =  PqT  is 
the  transmitted  bit  energy.  We  utilize  the  Gaussian  approximation  technique,  namely  we 
approximate  Z,  by  a  zero-mean  Gaussian  random  variable  with  variance  cr|  ,  which  denotes 
the  variance  of  Z,.  This  is  the  coherent  analog  to  the  Gaussian  approximation  found  in 
[5]  and  [6].  Given  this  approximation,  the  conditional  symbol  (bit)  error  probability  is 


Pc.i  =  2 


(18) 


in  which  erfc(  )  is  the  usual  complementary  error  function.  From  our  earlier  expressions, 
it  can  be  shown  that 

1-1 


2  ^ 

X 

Z,  2 


Eb 


+  9{\0Mk\)al, 

k:fil 


(19) 


where  is  the  variance  of  /fc,,.  If  we  once  again  assume  that  all  hits  that  occur  are  full 
hits,  we  find  that 

Eb 


1  ^ 
-  - 


2  UoJ 


k^x 


(20) 
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where  cr^  is  the  variance  of  the  interference  terms  under  the  assumption  that  any  hit  that 
occurs  is  a  fuU  hit.  It  can  be  shown,  using  the  results  of  [5],  that  aj  =  m,/,,  where  m,i.  is 
as  discussed  above. 

At  this  point  we  use  a  three-term  exponential  approximation  to  to  obtain 


3.1.3.  Packet  Throughput 

As  mentioned  earlier,  in  order  to  compute  packet  throughput  we  need  a  probability  of 
packet  correctness  or  erroi.  This  is  accomplished  by  deriving  first  the  conditional  packet 


error  probability  given  the  number  of  interfering  users  and  f ;  and  then  taking  expectations 
with  respect  to  the  distribution  of  interfering  (causing  hits)  users,  the  distribution  of 
potential  interfering  u.sers  in  the  cell,  and  the  distance  f  of  the  receiver  under  consideration 
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from  the  center  of  the  ceU.  The  first  step  is  straightforward;  we  can  obtain  the  desired 
conditional  probability  expressions  by  taking  the  known  error  probability  expressions  for  a 
noncoherent  FSK  or  a  coherent  BPSK  system  in  AWGN  from  the  literature  (e.g.  [9])  and 
substitute  for  the  variance  of  the  other-user  interference  and  the  three-term  exponential 
approximation  from  the  analysis  described  previously  in  this  section.  The  second  step  is 
demonstrated  below  for  the  cases  of  uncoded  packets  and  RS(n,it)-coded  (encoded  with  a 
Reed-Solomon  code  of  length  n  and  rate  k/n)  packets;  the  development  is  the  same  for 
other  coding  schemes. 

3. 1.3.1.  Uncoded  Systems 

For  a  general  uncoded  packet  of  length  L  we  have 

Pe.pkt  =  1  ~  • 

Since  Pc{f,K)  in  (11)  is  dependent  on  f,  on  the  number  of  potential  interferers,  and  on 
the  number  of  hits  given  A',  we  must  take  expectations  over  all  of  these.  Thus, 

P.,pkt  =  l-EfEKEK{PHf,P')}-  (23) 


where  q  is  the  number  of  different  frequencies  and  N,  the  number  of  symbols  per  hop. 
As  discussed  in  [7],  for  moderately  large  values  of  q,  (24)  either  bounds  or  approximates 


the  hit  probabilities  for  memoryless  random  hopping  patterns  and  Reed-Solomon  periodic 
hopping  patterns.  Thus,  (24)  is  a  good  general  expression  for  most  purposes. 

Suppose  now  that  the  summation  in  the  exponent  in  (11)  can  be  decomposed  into 
collections  of  interferers  having  equal  5(|0Af|).  We  would  then  have 

E  5(|0M*|)  =  E  (25) 

where  kj  is  the  number  of  interferers  with  p(lOA/|)  =  gj. 

Let  us  consider  some  planar  region  A  and  a  partition  of  A  given  by  such 


that  for  all  mobile  users  located  in  Aj,  <7(|0A/|)  =  gj.  Let  N{Aj)  denote  the  number 
of  such  mobile  users.  We  now  have  a  total  number  of  users  with  a  given  attenuation  and 
a  probability  that,  out  of  a  given  population  of  K  potential  interferers,  there  are  k  hits. 
Then  one  may  take  Kj  =  iV(Aj)  to  be  the  population  in  that  probability  expression. 
From  (11)  we  get  by  expanding  the  trinomial 


Pt{f,K) 


L  m 

EE 

m=0  t=0 


'L 


ti 


[Ci^(i,f)r[C2^(2,f)r 


X  [CsHid,  f)]^  exp  [[m^i  (f)  +  i^2(r) 


+(I_m-t)^3(r)]Ei?(|0^it|) 

k^i 


(26) 


Considering  this,  the  discussion  of  £/,{•},  and  the  fact  that  C,(f)  =  CiH{i,  f),  we  rewrite 


(23)  as 


(  \ 

/  \ 

L  m 

p..pkt  =  1  -  E  E 

L 

m 

m=:0  t=0 

\  '  > 

xEfEf^Eh 

{cr(r)C^(f)C 

X  exp([m^i(r)  +  J(52(^)  +  {L  -  m  -  i)^3{r)] 

xE.a)) 
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L  m 

=  i-EE 

m=0  t=0 


xEf{CT{f)C'^{r)C^-^-\f) 


X  exp  ([mii(f)  +  if)2{f)  ->r{L  -  m-  t)^3(r)] 


L  m 

=  i-EE 

m=0  1=0 

xf;,{c'r(f)c‘(f)c|'— ‘(f) 

^[•m6x(f)+iS2(f)+{L-m-i)S3(f)]g, 

U=1 

L  m 

=  i-EE 

tn=:0  t=0 

xEk 

Is*'  In  [Ph  exp  ([Tn6i(f )  +  i(52(f) 
+(i-m-i)63(r)]ff;)+ 1 -P/.]}}.  (27) 


Efc{-}  is  taken  by  following  a  technique  used  in  [1].  For  any  partition  of  the  surface  of 
the  cell  A,  say  defined  as  described  in  the  paragraph  following  (25),  and  constants 

{cj}"-, ,  we  have 


exp  ^A'(/l_,)aj  i  =  exp  A  5(/lj)  (c“^  -  1) 


where  5(>lj)  denotes  the  area  of  Aj.  Because  A'j  =  N(Aj),  (27)  and  (28)  can  be  combined 
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to  yield 


e,pkt  — 


/  \ 

(  \ 

L 

m 

K  y 

[  '  y 

'-EE 

tn=0  »=0 

xE,{c^{f)Ci{f)C^—'if) 


j=l 


X  exp 

+(L-m-i)63(f)]p_,)-  1])]}. 


(29) 


We  now  consider  the  limit  of  the  sum  inside  the  exponential  in  (29)  as  the  partition 
becomes  finer  and  finer;  it  approaches  an  integral  over  the  surface  A — standard 
arguments  for  Riemann  integration  can  be  used  here.  This  allows  Pe,pkt  to  be  written  as 


P^rki  =  1-EE 

m=0  »=0 


/  \ 

/  \ 

L 

m 

m 

\  / 

^  '  y) 

X  exp  l^L  (Ph  [exp([m^i(f)  +  i^2(r) 

+{L~m-  i)63if)]9{\0M\))  -  l])dHl]  } 


(30) 


where  denotes  the  integral  over  the  surface  A.  This  completes  the  analysis 

because  Ef{  /  is  simply  an  expectation  over  the  entire  area  A  under  consideration. 

3. 1.3. 2.  Coded  Systems 

We  now  consider  RS(Ti,A:)-coded  packets  with  error-only  correction  decoding.  The 
applicable  error  probability  expression  is  given  by 

(  \ 


Pe,RS  =  1-12 
t=Q 


P^-\\  -  Pc)' 


(31) 
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where  t  =  [n  —  k/2\  is  the  error  correction  capability  of  the  code.  This  can  be  also 
expressed  as 

With  respect  to  the  required  expectations,  this  can  be  treated  exactly  like  the  uncoded 
packet  error  probability  expression.  Therefore,  we  obtain 


t  t  n— Z+m  I  fi 

=  1-ee  e  e 

^=0  m=0  1=0  j=0  £ 


n  —  £  -h  m 


It  should  be  noted  that  in  the  case  of  coherent  BPSK,  the  three-term  exponential  approx¬ 
imation  must  approximate  P”,  rather  than  Pc,i,  when  considering  Reed-Solomon  codes. 
The  m  noted  here  refers  to  the  number  of  bits  in  BPSK  per  Reed-Solomon  symbol;  for 
example,  if  RS(32,16)  code  is  used,  m  =  5.  It  should  also  be  noted  that  for  the  Reed- 
Solomon  coded  systems  the  signal-to-noise  ratio  of  the  channel  symbols  is  that  of  the 
uncoded  system  multiplied  by  the  code  rate;  this  is  valid  under  the  assumption  that  the 
data  rate  of  the  system  remains  fixed  and  the  toted  bandwidth  of  the  system  expands  by 
the  inverse  of  the  code  rate. 

Once  one  has  obtained  an  expression  for  packet  error  probability,  the  expected  through¬ 
put  can  be  obtained  by  subtracting  it  from  one  and  multiplying  by  the  offned  traffic  load, 
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which  is  given  by  Gtot  =  A7rS(>l),  5(^4)  being  the  area  of  the  region  of  concern.  In  the 
cases  presented  here,  a  bandwidth  expansion  factor  must  also  be  introduced. 

'S  =  Gtot{l-Pc)B.  (34) 

B,  the  bandwidth  expansion  factor,  is  given  by 

B  =  -Rc,  (35) 

9 

where  q  is  the  number  of  hop  frequencies  and  the  code  rate  Rc  \s  1,  for  uncoded  packets, 
and  k/n,  for  RS{n,k)  coded  packets. 

3.2.  Bounds 

Error  probabilities  have  been  upper  bounded  giving  lower  bounds  on  system  through¬ 
put.  These  will  be  presented  here  for  the  purpose  of  comparison.  The  probability  of  a 
symbol  error  at  the  receiver  Pa(A'),  when  K  users  transmit,  can  be  upper-bounded  as 

p,{K)  <  1  -  (1  -  Po)(l  -  Phf-^  (36) 

where  Po  is  the  error  probability  for  a  single  user  system  in  the  presence  of  AWGN  with 
or  without  fading,  and  Ph  is  the  probability  of  a  hit  from  a  single  other  user;  expressions 
for  Ph  can  be  found  in  [7].  For  random  frequency-hopping  patterns,  we  can  use  Ph  = 
(1  -f  log2  A//A''fc)  where  Nb  denotes  the  number  of  bits  per  dwell  time  and  q  is  the 
number  of  frequencies. 

Let  us  discuss  now  the  calculation  of  throughput  when  the  total  number  of  users 
(packets)  transmitting  in  a  slot  (a  packet  slot,  which  is  L  symbols  long  in  the  uncoded 
case)  is  K .  Let  s{K)  and  Pe{K)  denote  the  average  number  of  packets  received  successfully 
and  the  average  probability  of  failure  (erroneous  reception)  of  a  typical  packet,  respectively, 
when  K  packets  are  transmitted.  This  gives  us  s{K)  —  K  -  KP^IK). 
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Let  A  be  Poisson-distributed  on  an  area  A  with  average  intensity  A  (packets  per  unit 
area).  Let  ^e(G)  and  P e{G)  denote  the  averages  of  the  quantities  s(A^)  and  Pg{K),  defined 
in  the  previous  paragraph,  with  respect  to  the  Poisson  distribution  p{K)  =  /  A'ij , 

where  i;{A^}  =  AA  =  G  (the  offered  load).  Then:P,(G)  =  l-[s,(G)/G]  =  [E{K Pg(K)} /G] 
and  the  normalized  throughput  is  given  by  t]{G)  =  G/q  [l  -  Pe(G)]. 

We  can  actually  show  that  P e(G)  =  E  {Pe{K  +  1)}.  Since  the  packet  error  probability 
Pe(A')  is,  through  the  symbol  error  probability  P,(A'),  a  function  of  A'  -  1  (the  number  of 
interfering  packets),  Pg(G)  is  a  function  of  K  and  not  of  K  —  1.  This  implies  that  in  the 
calculation  of  the  throughput,  instead  of  working  with  the  total  number  of  transmitted 
packets,  we  can  equivalently  work  with  the  number  of  interfering  packets  and  treat  that 
as  being  Poisson-distributed  with  the  same  parameter  G.  This  argument  also  justifies  the 
assumption  that  the  number  of  interfering  users  is  Poisson-distributed  with  parameter  G. 

We  can  now  write 

Pc(G)  =  1-P*(G)  =  P{[l-p,(A'+l)]}^ 

=  (1  -  Po)^e~'^(*-(’-^'‘>'']  (37) 

where  L  is  the  number  of  A/-ary  symbols  or  bits  transmitted  in  a  packet.  L  =  1  can  be 
used  to  bound  the  bit  error  probability.  The  throughput  can  be  then  lower-bounded  as 


V{G)  >  —(1  -  (38) 

q 

If  the  system  employs  error-control  coding,  we  can  follow  the  analysis  in  [11].  The 
analysis  will  be  shown  for  a  Reed-Solomon  (RS)  errors-only  correction  scheme;  other 
schemes,  such  as  RS  with  side  information,  have  been  analzved  similarly. 
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The  probability  of  a  symbol  error  Pa(A')  is  upper-bounded  in  a  way  similar  to  that 
in  (36),  where  P/,  is  the  probability  of  a  hit  for  an  RS  symbol,  consisting  of  me  m-ary 
symbols;  thus,  we  must  use  me  log2  A/  instead  of  log^  M  in  the  expression  that  provides 
the  probability  of  a  hit.  When  the  RS{n,k)  code  errors-only  correction  is  employed  and 
one  codeword  per  packet  is  transmitted,  the  probability  of  correct  reception  of  a  codeword 
is  given  by 

/  \ 


Pc{K)  > 


[p,(A')?  [i-p,(A')r^ 


(39) 


where  i  =  [n  —  k/2\  is  the  error-correction  capability  of  the  code  and  P,(A  )  <  1  - 
(1  —  Finally,  we  can  get  an  expression  for  the  average  probability  of  correct 

packet  reception  Pc{G)  that  allows  us  to  wTite  the  following  expression  for  the  normalized 
throughput  (per  frequency  slot)  t/(G) 

viG)  > 

X  [(l-Po)'"^(l-P/.)'^'J"'' 


\ 

'  y 

4.  ACCURATE  EVALUATION  OF  BIT  ERROR  PROBABILITY 
FOR  COHERENT  MODULATION 


Gky^^ 
q  n 

”  J=:0  «=0 


1 


(10) 


If  coherent  modulation  is  employed,  the  resulting  system  can  be  analyzed  by  means  of 
a  characteristic  function  method,  as  will  be  demonstrated  below. 

In  this  case,  the  system  will  be  assumed  to  utilize  binary  pha,se-shift  keying  (BPSK) 
with  coherent  demodulation.  The  signal  transmitted  by  the  Arth  mobile  station  is  of  the 
form  given  in  (12).  We  assume  here  that  Rayleigh  fading  may  or  may  not  be  present.  We 
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modify  (13)  to  get 


r*(0  =  F,^2P{M,)h[t)^[t)  cos[2:r[/,  +  +  h{i)]  (41) 

where  Fk  is  a  Rayleigh  fading  factor,  if  fading  is  assumed  to  be  present  and  all  other 
quantities  are  as  previously  described. 

The  dehopper  output  rj(<)  used  in  (15)  is  now  given  by 

U(t)  =  '^\^Fksl^^nfk{t-rk)JM)]bk{t-rk)^{t-Tk) 

X  cos[2TTfat  +  (t>k(t)]} Vd(t]  (42) 

We  proceed  until  we  again  have  a  modified  test  statistic  Z,  of  the  form  given  in  (17).  Our 
goal  is  to  derive  expressions  for  the  probability  of  error  and  the  expected  throughput  q. 
Since  a  coherent  receiver  is  being  used  here,  the  characteristic  function  method,  described 
in  [13]  and  simplified  to  the  case  where  the  signature  sequences  are  identically  1.  will  be 
used  to  derive  the  desired  probability  expression. 

The  computation  of  the  error  probability  involves  the  characteristic  function  of  the 
noise  h,  which,  in  our  case,  is  $n(«)  =  exp  (  — u^/2  •  A’o/2£'i,,o).  and  the  characteristic 
function  of  the  interference  process,  which  is  the  middle  term  in  (17).  This  will  be  derived 
here  for  the  synchronous  case  and  for  the  asynchronous  case,  both  with  and  without 
Rayleigh  fading. 

4.1.  Synchronous  Systems 

In  the  synchronous  frequency-hopping  ca.sc.  we  have  =  0  for  all  k.  The  kth  signal 
interferes  with  the  i-th  signal  if  and  only  if  the  two  signals  are  exactly  superimposed  upon 
each  other;  this  suggests  that  user  k  interferes  with  user  i,  if  there  is  a  full  hit  when  two 
signals  are  sent  in  the  same  frequency  “bin”  at  the  same  time.  Thus,  if  q  denotes  the 
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number  of  different  frequency  bins,  the  probability  of  a  hit  is  I/9. 

A  result  used  in  [l]  can  now  be  adapted  to  give  an  expression  for  ^i{u),  the  character¬ 
istic  function  of  the  interference  process.  The  result  used  in  [1]  is  that,  if  the  interfcrf'nce 
process  can  be  written  as  /  =  for  ^  random  process  {A^}  with  pdf  pxv') 

and,  if  the  number  of  users  is  Poisson  distributed  on  the  area  A  with  parameter  (mean 
intensity)  A,  then 


^liu)  =  exp 


l>/.{|/ 


px{x)e"‘'^^^^dx 


-  1 


(43) 


in  which  /(A/)  M  denotes  a  surface  integral  and  A  represents  the  area  of  the  cell 
(area  of  concern).  This,  however,  is  not  the  exact  form  that  must  be  used  here,  since 
the  interference  process  in  our  problem  has  the  form  of  the  middle  term  of  (17)  and  the 
probability  of  hits  must  also  be  taken  into  account.  As  a  result,  the  characteristic  function 
of  the  interference  process  takes  the  form 


$/(u) 


dHl 


(44) 


where 


1  7^" 

4>(u)  =  —  /  exp{iucos9)d0  =  do(u).  (45) 

ZTT  Jo 

In  (45),  the  cos  0  term  was  obtained  by  considering  (43)  and  (17),  letting  r/;  =  0  and 
assuming  uniformly  distributed  random  phase  processes  {6  is  uniform  on  the  interval 
[0,27r));  Jo{-)  is  the  Bessel  function  of  the  first  kind  and  zeroth  order. 


If  we  consider  our  region  of  concern  to  be  a  disk  of  radius  R,  then  the  integral  in 
(44)  can  be  easily  computed.  Let  us  denote  by  r  the  radial  distance  |OA/|  and  con¬ 
sider  the  annulus  (r,r  -f  dr).  Then  we  can  write  that  [do  (uv/sf(|OA/|)j  -  ij  d^M  = 
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2xr  (^\/30l)  -  l]  dr,  where  ff(-)  is  as  in  (17).  This  expression  can  be  evaluated 
using  formulas  found  in  [14]  to  yield  series  solution. 

Let  us  consider  the  same  model  with  Rayleigh  fading.  This  is  accounted  for  by  means 
of  Fk  in  the  expressions  above.  The  fit  are  assumed  to  be  independent  and  identically 
distributed  with  a  Rayleigh  parameter  6^;  therefore,  F{F^}  =  2b^.  As  a  result  of  the 
addition  of  fading,  (45)  becomes 


yoo  I  y2)r 

$(u)  =  /  iV(y)r-  /  exY>{iuy  cos  e)d9dy 

Jo  Jo 


=  exp( - ~) 


(46) 


In  (46),  Pfiy)  —  exp  (— y^/26^)  is  the  Rayleigh  density  function  with  parameter 

b^.  Since  the  product  of  a  Rayleigh  random  variable  and  the  sine  or  cosine  of  a  uniform 
on  [0,27r)  random  variable  yields  a  Gaussian  random  variable,  G  =  Fcos9  has  pdf 
Pair])  =  l/(\/27r6)  exp  (-77^/26^)  and  the  right-hand  equality  is  derived  by  performing 
the  integration  with  respect  to  p. 

If  the  circular  region  model  mentioned  above  is  used,  (44)  becomes 


{ X  /•«  / 

$/(u)  =  exp  {—I  j  e: 

[q  Jo  \ 

=  exp|^([e-«<“)^o''_i 

,  r  [e--'-  - . 


b’^u'^ 


—  1  1  2'KrdT 


■ktI 


2-Krdr 


(47) 


where  B{u)  =  b^u'^/2.  A  closed-form  result  can  be  obtained  using  a  formula  from  [16], 
if  o  is  known.  If  the  typical  value  a  =  3.5  is  used,  (47)  becomes 


^l(u) 


exp  (7iT^[c  ^*“>’■0  -  1]  -  7r[/?^  -  rl] 
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I  1 ^  ~  ^  T-6-it 


(48) 


4.2.  Asynchronous  Systems 

We  now  focus  on  the  asynchronous  frequency-hopping  case.  It  will  be  assumed  that 
the  Tk  in  (42)  are  independent  and  identically  uniformly  distributed  on  the  interval  [0,  T/,]. 
As  a  result  of  these  delays,  there  are  three  possible  interference  scenarios  for  a  signal 
transmitted  in  a  given  frequency  bin:  no  hit,  partial  hit,  or  full  hit.  These  must  all  be 
accounted  for  in  the  computation  of  the  characteristic  function  of  the  interference  process. 

Since  the  interference  term  in  (17)  is  a  sum  of  independent  random  variables,  the 
overall  characteristic  function  can  be  seen  as  a  product  of  the  individual  characteristic 
functions,  that  is. 


4*/(ti)  =  ElU 


j 


=  f;|n[i-F  -b  P/$^*’’'ku)  +  Pp$(*'')(u)]| 


(49) 


In  (49),  E  denotes  expectation  with  respect  to  the  number  of  interfering  users,  Pf  is 
the  probability  of  a  full  hit,  Pp  is  the  probability  of  a  partial  hit,  1-  P  =  I  -  Pj  —  Pp  is  the 
probability  of  no  hit  occurring,  and  $y*’’^(w)  and  $p*’'*(u)  are  the  characteristic  functions 
of  Ik,i,  given  that  a  full  hit  or  a  partial  hit  hais  occurred.  Note  that  for  the  case  of  no  hit, 
h,t  =  0,  which  results  in  a  characteristic  function  equal  to  1.  The  validity  of  the  second 
equality  in  (49)  was  proven  in  [13]. 

The  quantity  inside  the  expectation  in  (49)  can  be  written  in  the  form 


expjEM.ln  [l  -  ^  +  +  Pp#‘’(«)]}. 


Let  us  consider  a  step  function  g{\OM\)  =  E"=j  9}  where  is  a  partition 

of  the  region  A;  /s(M)  is  1,  if  M  G  B,  and  0  otherwise;  and  g{\OM\)  =  for  all 


24 


M  g  Ay  If  we  denote  by  N{Aj)  the  number  of  users  in  the  area  Ay  we  obtain  that 


$/(«)  =  E  |exp  A^(^,)ln  [l  -  P  +  P/$/(u^)  + 


(50) 


in  which  $/(uy/5J)  and  denote  the  characteristic  functions  of  the  interference 

due  to  any  particular  user  located  in  for  the  cases  of  full  hits  and  partial  hits,  respec¬ 
tively.  Then  the  expectation  can  be  written  as 


$/(u)  =  exp|APyX:5(yl,)[$/(u^)-l] 

+  APpf^5(>l,)[$p(uv^)-l] 
j=i 

In  the  general  case,  g{\OM\)  is  a  limit  of  step  functions,  implying  that  $/(uv/p(l0M|)) 
and  ^p{u\/g{JOMj))  are  such  limits,  too,  and  the  overall  characteristic  function  becomes 


$/(u)  =  exp^AP/  j  -  1 

+  ^PvJ^  [^P  (u^pdOMI))  -  l]  d^Mj  . 


(52) 


Expressions  for  Pj  and  Pp  were  given  in  [13]  for  first-order  Markov  random  hopping 
patterns,  for  memoryless  random  hopping  patterns,  and  for  Reed-Solomon  periodic  hop¬ 
ping  patterns.  Following  [13]  for  the  case  in  which  the  signature  sequence  is  identically  1, 
we  obtain  that 
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1  dr 
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in  the  absence  of  fading,  and 
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in  the  presence  of  Rayleigh  fading.  In  (53)  and  (54),  the  functions  and  R.^  are  the 
partial  autocorrelation  functions  of  the  shaping  waveform  ^’(0  and  are  defined  by  R^{s)  = 
-  s)dt  and  R^{s)  =  R^{T  -  s). 

The  relationship  in  [13],  which  expresses  the  bit  error  probability  of  the  receiver  as 
an  integral  involving  the  characteristic  functions  of  the  noise  and  the  inverference  process 
(see  [6]),  can  be  modified  to  depend  on  the  position  the  “transmitter  of  concern”  in  the 
following  way 


n(lOMl)  =  Q{J2Eb,o/No 


sin[^(|(>M))w] 


$fi(u)  [1  -  ^►/(t/)]  du. 


where  |0M|  represents  the  distance  between  the  receiving  (base)  station  and  the  “trans¬ 
mitter  of  concern”  (the  ith  user  in  the  analysis  above),  Q{x)  is  the  complementary  error 
function  given  by  Q(x)  =  exp  (-y^/2)  dy,  and  $fi(u)  was  defined  earlier.  $/(u) 

is  given  by  (44)  for  synchronous  systems,  by  (45)  for  no  fading,  and  (46)  for  Rayleigh  fad¬ 
ing.  For  asynchronous  systems,  it  is  given  by  (52)  and  (53)  for  no  fading  and  by  (54)  for 
Rayleigh  fading. 

The  result  in  (55)  is  useful  because  it  gives  us  a  measure  of  the  correctness  of  each 
individual  bit.  Unfortunately,  to  evaluate  the  throughput  one  needs  the  packet  error 
probability  Pe(|OA/|),  which  can  not  be  expressed  as  an  explicit  function  of  Pt,{\OM\). 
The  error  probability  of  each  one  of  the  L  bits  in  a  packet  depends  on  the  number  of 
interfering  users  and  their  received  power  at  the  base  station.  Therefore,  one  should  first 
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derive  the  bit  error  probability  given  the  number  of  interfering  users;  from  that  one  should 
derive  the  packet  error  probability,  given  again  the  number  of  interfering  users,  and  then 
finally  average  the  latter  with  respect  to  the  Poisson  distributed  number  of  interfering 
users.  As  a  restilt,  the  best  that  can  be  done  using  this  method  for  packets  (coded  or 
uncoded  ones)  is  an  approximation. 

5.  NUMERICAL  RESULTS 

In  discussing  our  numerical  results,  we  must  note  several  things.  First,  we  take  A 
to  be  a  circular  region  of  radius  Ji  and  then  we  set  R  =  1.  This  can  be  viewed  as  a 
normalization  and  the  results  can  be  generalized  to  regions  of  larger  or  smaller  size;  this 
is  especially  evident  if  one  considers  that  the  expected  number  of  users  in  A  is  given  by 
N  =  XitR?  and  that  our  results  are  indexed  by  N  (and  therefore  by  A). 

A  brief  discussion  of  the  issue  of  normalizing  the  radius  R  of  the  cell  is  in  order  here. 
We  assume  that  in  the  mobile-to-base  mode  of  operation,  all  mobiles  transmit  with  the 
same  signal-to-noise  ratio  Ej/Aq,  which  is  such  that  the  signal-to-noise  ratio  received  at 
the  base  from  a  mobile  at  a  distance  r,  given  by  =  (Efc/Ao)  ff(r),  is  at  least  as  large 

as  EbJNo  =  {Eb/No)g{R),  the  received  energj'  from  a  mobile  at  the  periphery  of  the  cell. 
Therefore,  if  Eb/No  is  set  at  the  minimum  acceptable  level  for  reliable  communications, 
the  transmitted  signal-to-noise  ratio  for  all  mobiles  is  given  by  Eb/No  =  ^Ffc/Ao)  /ff(iZ); 
if  i?  =  1,  Eb/No  =  Eb/No,  while  if  iZ  >  1,  g{R)  can  be  computed  and  used  to  compute 
Eb/No  from  the  above  formula,  while  using  /Z  =  1  for  all  other  calculations. 

It  is  also  important  to  consider  the  computational  requirements  of  our  results.  Looking 
at  (30)  or  (33),  one  can  see  that  there  are  multiple  levels  of  summations  and  integrals 
involved  in  the  computation  of  each  value  of  the  packet  error  probability  (note  that  Ff  {  •}  is 
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an  integral).  In  computing  the  results  shown  here,  we  sought  to  avoid  this  last  expectation 
over  f  and  explored  two  alternative  methods. 

The  first  method  was  to  set  f  equal  to  some  value  such  that  a  given  probability  of 
acceptable  communication  is  maintained.  As  discussed  in  [1],  it  can  be  shown  that  in 
this  case,  r  =  y/PR,  where  P  is  the  probability  of  acceptable  communication  to  be 
liiaintaired.  If  one  considp»‘c  vabies  of  P  clo,ce  to  one — a-s  would  be  tb“  range  of  values  one 
might  like  to  examine — one  can  see  that  the  results  obtained  will  be  overly  pessimistic. 

As  a  result  of  the  pessimistic  nature  of  the  first  method,  we  have  opted  to  use  another 
method,  namely  to  set  f  equal  to  a  value  such  that  an  average  value  of  g(f)  is  obtained; 
then  f  =  where  a  and  g(-)  are  as  in  (2).  We  find  that  the  expectation  is 

given  by 

{i/C'-)}  =  ^  <“) 

The  value  of  this  method  is  due  to  the  fact  that  the  probability  of  symbol  or  packet 
correctness  is  usually  a  convex  n  function  of  g(r)  in  part  of  the  range  of  r.  In  this  range, 
we  can  use  Jensen’s  inequality 

^r{Pc(ff(r))}  <  Pc{Er(ff(r))}  =  Pc(g(f))  (57) 

to  obtain  the  \ipper  bound  on  the  right-hand  side  which  depends  only  on  g(r).  In  the 
range  of  r  in  which  the  function  Pc  is  not  convex  n  in  r,  we  can  still  use  Pc(g(f)  as  an 
approximation. 

In  all  the  numerical  results  presented  in  this  section  the  signal-to-noise  ratio  (denoted 
by  SNR  in  the  figures)  is  the  Eb/No  with  respect  to  the  AWGN  defined  in  the  second  para¬ 
graph  of  this  section.  The  number  of  bits  per  hop  (dwell-time)  is  Nb  =  10,  which  implies 
that  in  the  M-ary  FSK  case  there  are  two  32-ary  symbols  (of  5  bits  each)  transmitted  in 


each  hop.  Moreover,  each  symbol  of  the  RS  codes  with  codeword  lenght  n  =  32  contains 
5  bits  and  thus  two  RS  symbols  are  transmitted  on  one  hop.  The  number  of  frequencies 
available  for  hopping  is  q  =  100. 

Fig.  2  compares  the  results  using  the  Gaussian  approximation  technique,  described  in 
Section  III-A,  with  the  bound  presented  in  Section  III-B  for  the  case  of  coherent  signaling 
and  RS(32,  16)  coding  v'ith  error-only  correction  decoding  This  demonstrates  that  the 
bound  tends  to  be  extremely  pessimistic.  On  the  other  hand,  the  Gaussian  approximation 
tends  to  be  moderately  optimistic  for  most  of  the  range  of  N  (the  average  number  of 
neighbors).  This  last  argument  is  justified  as  follows;  Gaussian  approximations  tend  to  be 
usually  optimistic  when  applied  to  multiple-access  problems  (refer  to  [5]  and  [6]).  However, 
the  assumption  that  aU  interfering  users  cause  full  hits  is  pessimistic.  Furthermore,  the 
use  of  g(f)  results  in  pessimistic  values  of  the  probability  of  correct  reception  for  part  of 
the  range  of  r.  Thus,  overall  the  Gaussian  approximation  used  here  is  expected  to  be 
moderately  optimistic.  This  can  be  actually  verified  through  a  comparison  to  the  exact 
expressions  in  Fig.  6  and  7. 

Fig.  3  deals  with  the  coherent  case  and  demonstrates  the  decrease  in  throughput  with 
decreasing  Reed-Solomon  code  rate.  As  one  might  expect,  throughput  is  proportional 
to  code  rate  due  to  the  normalization  by  the  code  rate.  However,  the  RS(32,8)  code 
provides  slightly  more  than  half  of  the  RS(32,16)  throughput  when  N  =  150,  which  is 
explained  by  an  improvement  in  the  expected  probability  of  error.  In  both  Fig.  2  and  3, 
the  values  of  the  throughput  for  large  N  (larger  than  150)  are  optimistic  and  should  not 
be  considered  fully  dependable.  This  is  due  to  the  fact  that  the  accuracy  of  the  Gaussian 
approximation  used  in  the  calculation  of  the  conditional  probability  of  error  in  Section 


29 


3.1.1  degrades  gradually  a.s  the  number  of  interfering  users  increase;  actually,  it  becomes 
overly  optimistic  for  a  large  number  of  interfering  users. 

We  now  move  on  to  the  noncoherent  case  (32-ary  FSK  with  noncoherent  demodula¬ 
tion).  Fig.  4  corresponds  to  Fig.  2.  Once  again  it  is  demonstrated  that  the  bound  is  very 
pessimistic  and  the  Gaussian  approximation  overly  optimistic  for  large  N . 

Fig  5  presents  results  using  RS(32,  16)  and  RS(32,  8)  coding.  The  same  observations 
that  were  made  regarding  Fig.  3  should  be  made  here. 

Fig.  6  demonstrates  the  relationships  among  the  exact  expression  developed  in  this 
section  and  the  Gaussian  approximation  and  bound  presented  in  Section  III.  The  bound 
tends  to  be  very  pessimistic.  The  Gaussian  approximation,  on  the  other  hand,  tends  to  be 
somewhat  pessimistic  for  low  to  moderate  values  of  N  (the  average  number  of  neighbors), 
while  it  provides  a  closer  approximation  and  even  becomes  slightly  optimistic  for  moderate 
to  K  -ge  values  of  N.  This  reflects  on  the  previous  throughput  results:  for  small  TV,  the 
approxim  tion  is  closer  to  the  exact  value  than  the  bound  but  is  pessimistic;  for  moderate 
TV,  it  is  close  ’.o  the  exact  value;  for  large  TV,  the  results  exceed  what  they  should  be  and 
this  effect  is  magmt'pd  by  the  RS  coding. 

Fig.  7  shows  results  for  synchronous  and  asynchronous  systems.  These  results  are  so 
close,  that  they  appear  as  <  single  curve  in  the  figure.  The  closeness  of  the  results  for 
the  synchronous  and  asynchronous  cases  is  due  to  the  fact  that  the  probability  of  a  hit 
IS  Ph  =  l/q  for  synchronous  systems  and  P/,  =  (1  -|-  1/TV(,)  /q  for  asynchronous  systems, 
whereas  m^,  =  1/3  for  asynchronous  systems  and  =  1/2  for  synchronous  systems, 
thus  balancing  out  the  effect  of  other-user  interference  of  the  two  systems.  Also  shown  is 
a  curve  that  represents  an  asynchronous  system  with  Rayleigh  fading;  the  parameter  of 
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Rayleigh  fading  is  =  0.01,  which  implies  that  the  desired  signal's  SNR  is  attenuated  by 
.02  on  the  average,  that  is  E{SNRjad}  =  26^  ■  SNR  =  .02  •  SNR.  It  is  assumed  that  the 
Rayleigh  fading  is  sufficiently  slow  to  make  coherent  demodulation  feasible.  As  expected, 
this  system  does  significantly  worse  than  its  counterpart  without  fading. 

6.  CONCLUDING  REMARKS 

We  presented  methods  by  which  error  probability  and/or  throughput  can  be  computed 
for  FH  cellular  mobile  packet  redio  networks.  These  methods  encompass  a  Gaussian  ap¬ 
proximation  applicable  to  all  systems  and  an  exact  method  that  can  be  applied  only 
to  coherent  systems.  It  was  shown  that  the  previously  derived  lower  bound  on  system 
throughput  (upper  bound  on  error  probability)  is  a  loose,  pessimistic  bound  in  compar¬ 
ison  to  the  approximation  and  the  exact  results.  The  approximation  was  shown  to  have 
satisfactory  accuracy  for  the  lower-to-moderate  values  of  the  average  number  of  neighbors, 
but  to  be  overly  optimistic  for  large  values  of  the  average  number  of  neighbors. 

One  major  problem  encountered  lay  in  the  computational  complexity  of  the  expres¬ 
sions.  In  general,  each  curve  in  the  figures  presented  required  a  considerable  amount  of 
CPU  time  to  compute.  In  the  case  of  some  of  the  exact  expressions  (especially  the  asyn¬ 
chronous  case  without  fading),  this  computation  was  extremely  long.  We  used  a  VAX  780 
computer  and  a  few  hours  of  supercomputer  CPU  time.  What  this  indicates  is  that  this 
work  may  very  well  represent  a  practical  limit  as  to  how  precisely  one  can  characterize 
these  systems  and  still  be  able  to  compute  analytically  the  various  performance  measures 
of  interest. 
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Throughput  versus  average  number  of  neighbors  for  varying  code  rate  (Reed-Solonion  coding 
with  error-only  correction  deccxiing,  coherent  BPSK,  SNR  =  15  dB,  q=l(X),  N  /i,  =  10 
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Fig.  4  —  Throughput  versus  average  number  of  neighbors:  Gaussian  approximation  versus  bound 
(RS  (32. 16)  ccxling  with  error-only  correction  decoding,  noncoherent  32-ary  FSK,  SNR=  15  dB,  q=  I(X)  N 


5  —  Throughput  versus  average  number  of  neighbors  for  varying  code  rate  (Reec-Solomon  ctxling 
with  error-only  correction  decoding,  noncoherent  32-ary  FSK,  SNR  =  15  dB,  q=  100,  N  *  =  10). 
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Fig.  6  —  Bit  error  probability  versus  average  number  of  neighbors:  exact  method,  Gaussian 
approximation,  and  bound  (coherent  BPSK,  asynchronous,  SNR=  10  dB,  q=  1(X),  N  h=  10). 


Fig.  7  —  Bit  error  probability  versus  average  number  of  neighbors;  synchronous  versus  asynchronous  md  faded 
versus  not  faded  (coherent  BPSK,  SNR=  10  dB,  q=  1(X),  N  j,  =  10,  b^  =0.01  when  fading  is  present). 


